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ABSTRACT 


An analysis is presented of the failed re-development of ex-Gaston and the development of 
Karl during the 2010 PREDICT field campaign based on the European Centre for Medium 
Range Weather Forecast (ECMWF) analyses. This study analyzes the flow kinematics and 
dynamics of ex-Gaston to investigate the role of dry, environmental air in the failed redevel¬ 
opment. The findings support the hypothesis that entrained, dry air near 600 hPa thwarted 
convective updraughts and vertical mass flux, which in turn leads to a reduction in vorticity 
and a compromised pouch at these middle levels. A compromised pouch allows further 
intrusion of dry air and inhibits vorticity amplification. This study supports recent work 
investigating the role of dry air in moist convection during tropical cyclogenesis. 

The circulation and vorticity dynamics of Karl are used to describe the system-scale spin- 
up within the marsupial paradigm framework. The circulation tendency shows that Karl 
underwent periods of spin-up and spin-down until its eventual development on 14 Septem¬ 
ber. These spin-up and spin-down episodes correspond with the cycle of convective activity 
in the Caribbean. The development process is concomitant with recent work describing the 
genesis and intensification process in developing tropical cyclones. 
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CHAPTER 1: 
Introduction 


A major, unsolved problem in dynamical meteorology has been the genesis of tropical cy¬ 
clones, hurricanes, and typhoons. A fundamental problem in understanding tropical cyclo¬ 
genesis is describing the interaction of the approximately conserved properties of potential 
vorticity and moist entropy. In the innovative work of Dunkerton et al. (2009), hereafter 
DMW09, the authors propose a new framework for understanding these interactions in 
tropical cyclogenesis. In their description of tropical cyclogenesis, DMW09 use a La- 
grangian viewpoint to describe the interaction of an easterly tropical wave, the synoptic 
environment, and deep convection. As the wave propagates primarily westwards, the flow 
establishes dynamic flow boundaries (a Kelvin cat’s eye) that effectively trap moist parcels. 
The region where this trapping occurs is referred to as the “wave-pouch." This pouch is cen¬ 
tered about the intersection of the wave trough-axis and the critical latitude (Figure 1.1). 
The critical latitude is the latitude where the mean flow and wave phase speeds are equal 
(DMW09). 

This new framework is referred to as the marsupial paradigm. The analogy being that a 
marsupial protects its young in a pouch, safe from a hostile environment. As an easterly 
wave propagates westward, it can provide a region where an incipient proto-vortex is pro¬ 
tected from hostile environmental conditions, such as vertical shear or dry air. The tropical 
wave supports the growth of the proto-vortex until it has intensified and is self-sustaining. 
In a study by Wang et al. (2009), the authors found that there was a “sweet spot" (Fig¬ 
ure 1.1) or preferred location for genesis. The sweet spot is the intersection of the critical 
latitude and the wave trough-axis described by DMW09. 


1 





Figure 1.1: Schematic depicting the wave-pouch region. The dashed-green lines are the 
horizontal wind currents of the easterly wave in a Eulerian reference frame. The black 
curves represent the Lagrangian wave-pouch boundary as revealed in the frame moving 
with the frame of the easterly wave disturbance. The purple line depicts the critical latitude 
(CL). The straight black line labeled “Trough" represents the wave-trough axis. The red 
arrow is an easterly jet. SAL is the dry, Sahara Air Layer. The arrow labeled “Inflow" 
shows a potential pathway into the wave-pouch (from Montgomery 2010a). 


Lramed by this marsupial paradigm, DMW09 pose three, testable hypotheses for tropical 
cyclogenesis: 


1. The cat’s eye region within the critical layer has cyclonic recirculation with 
little strain deformation, and its center serves as the focal point for vorticity 
aggregation of convectively generated cyclonic vertical vorticity and their vor¬ 
tical remnants. 

2. Air inside the cat’s eye is repeatedly moistened by convection and is pro¬ 
tected to some extent from the lateral intrusion of dry air. 

3. The parent wave is maintained and possibly enhanced by convectively am- 


2 






plified mesoscale eddies within the wave. (Rutherford and Montgomery 2012) 


The new cyclogenesis model and accompanying scientific hypotheses were established ob- 
servationally in the Atlantic and eastern Pacific sectors by DMW09. 

In order to get a better understanding of tropical cyclogenesis and test the marsupial 
paradigm, the Pre-Depression Investigation of Cloud Systems in the Tropics (PREDICT) 
experiment was conducted during the 2010 Atlantic hurricane season (Montgomery et al. 
2012). This experiment was designed to identify tropical waves before they were identified 
as tropical depressions. The environmental characteristics of these mesoscale systems were 
measured primarily by dropsonde delivered by a Gulfstream V (GV). The aircraft was op¬ 
erated from St. Croix, U.S. Virgin Islands, and had an operational range from the Atlantic 
Basin to Central America at a range of cruising altitudes between 12 km to 14 km. 

PREDICT Mission scientists used “pouch tracking" techniques (Boothe 2015a) based on 
the marsupial paradigm to identify potential systems. Pouch tracking involves using model 
forecasts to analyze the evolution of tropical disturbances from a Lagrangian perspec¬ 
tive, as well as environmental properties like vertical shear and relative humidity. This 
method allows forecasters to predict the likelihood of the development tropical distur¬ 
bances. The PREDICT experiment investigated eight tropical systems; two nondeveloping 
African waves, three developing waves, one developing Caribbean gyre, one failed redevel¬ 
opment of an African wave, and one system of undetermined fate (Montgomery et al. 2012, 
Table 2 therein). There were a total of 25 research flights over 24 days with 561 dropsondes 
deployed. The PREDICT experiment provided an unprecedented amount of in-situ obser¬ 
vations of developing and nondeveloping tropical systems. Concurrent with the PREDICT 
experiment, NASA and the National Oceanic and Atmospheric Administration (NOAA) 
were also conducting experiments to study the tropical disturbances. The NASA Genesis 
and Rapid Intensification Processes (GRIP) used a DC-8 aircraft to deploy dropsondes in 
tropical disturbances. The NOAA Intensity Forecasting Experiment (IFEX) was primarily 
focused on observing developed systems, and used a G-IV to deploy dropsondes and two 
WP-3 aircraft to deploy dropsondes and perform airborne Doppler radar observations. 

One of the first studies of PREDICT data was performed by Smith and Montgomery 
(2012). In their study, the authors characterized the convective environments of three sys- 
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terns. Two of the systems underwent genesis: pre-Karl developed into Hurricane Karl and 
pre-Matthew developed into tropical storm Matthew. The other system, named Gaston, 
achieved tropical storm status briefly, but eventually weakened and failed to redevelop. 
Because this system never redeveloped into a tropical storm, it is referred to hereafter as 
ex-Gaston. Smith and Montgomery (2012) calculated the convective available potential 
energy (CAPE), convective inhibition (CIN), and moisture properties. They found that the 
lower troposphere tended to become progressively more moist, but in the case of ex-Gaston 
the troposphere became drier between an altitude of 2 and 9 km. Based on these findings, 
Smith and Montgomery (2012) hypothesized that dry air (or “anti-fuel") weakened the con¬ 
vective updraughts of ex-Gaston and limited the amplification of absolute vorticity, thereby 
preventing the system from redeveloping. Sauer (2014) performed a complimentary study 
where he analyzed the dropsonde data for the PREDICT systems not reported in Smith and 
Montgomery (2012). Komaromi (2013) analyzed composite dropsonde data for the de¬ 
veloping and non-developing system during PREDICT and had similar findings as Smith 
and Montgomery (2012), namely progressive moistening during the developing cases and 
progressive drying during the non-developing cases. 

Wang (2012) provided additional evidence for a progressive moistening of the pouch in 
developing cases, and drying in non-developing cases. Wang (2012) used numerical sim¬ 
ulations of the development of Hurricane Felix (2007), studied earlier by Wang et al. 
(2010a,b), using the Weather Research and Forecasting (WRF) Model and PREDICT data 
to examine thermodynamic variables as a function of radius from the disturbance’s sweet- 
spot. She found that for developing tropical cyclones, the inner pouch region is moist 
compared to the outer pouch region, and becomes progressively more moist as the sys¬ 
tem organizes. Based on the results of Wang (2012), it seems reasonable to attribute the 
favorable thermodynamic conditions to the minimal shear and strain within a protective 
pouch. 

Davis and Ahijevych (2012) used the PREDICT data to study the structural evolution of ex- 
Gaston, Karl, and Matthew (the same systems analyzed by Smith and Montgomery 2012). 
They found that ex-Gaston started as a moist core surrounded by a dry environment, and 
vertical shear caused the circulation centers to be vertically misaligned. Davis and Ahi¬ 
jevych (2012) conclude that the dry air limited the convective activity by enhancing the 
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potential for downdraughts near ex-Gaston’s pouch. They also found that Karl’s circula¬ 
tion centers were initially vertically misaligned, but aligned later during its evolution. After 
the alignment, Karl underwent genesis within a moist environment. Matthew also formed 
within a moist environment, and was vertically aligned throughout its relatively rapid de¬ 
velopment. Davis and Ahijevych (2012) concluded that the marsupial paradigm helps to 
explain a level of predictability of tropical cyclogenesis. 

In a study by Rutherford and Montgomery (2012), the authors introduced a new method for 
analyzing tropical cyclogenesis. They examined the flow topology and flow boundaries in 
ex-Gaston and Karl using the ECMWF model data. The Lagrangian methods presented by 
Rutherford and Montgomery (2012) test the first and second hypotheses of the marsupial 
paradigm. First, the authors were able to calculate the flow boundaries and identify a 
pouch. Second, Rutherford and Montgomery (2012) used particle trajectories to identify 
possible dry air intrusion. Third, they showed evolution of relative vorticity within the 
pouch for ex-Gaston and Karl. Rutherford and Montgomery (2012) showed that ex-Gaston 
had a pouch that was open to hostile, dry air which they argued led to this system’s failed 
redevelopment. Further, they found similar results as Davis and Ahijevych (2012), namely 
Karl initially had a vertically misaligned pouch, that eventually aligned, and the system 
underwent genesis. 

The marsupial paradigm finds additional support in idealized numerical modeling simula¬ 
tions (Montgomery et al. 2010b; Nicholls and Montgomery 2013; Wang et al. 2010a,b), 
recent case studies in the field in the western North Pacific during the Tropical Cyclone 
Structure Experiment 2008 (TCS08, Fussier III 2010; Fussier III et al. 2014; Montgomery 
et al. 2010a, 2012; Raymond and Fopez-Carrillo 2011), in NASA’s ongoing Hurricane and 
Severe Storm (HS3) missions (2012-15), and the case of Hurricane Sandy (Fussier III et al. 
2015). The field data afford a resolved view of horizontal and vertical structure in the wave 
pouch and its immediate surroundings, valuable for system centering, circulation magni¬ 
tude, vorticity balance, interleaving of air masses, and moist thermodynamic profiles. 

A somewhat different perspective on tropical cyclogenesis is offered by Raymond and Ses¬ 
sions (2007) and Raymond and Fopez-Carrillo (2011). In a review of the numerical ex¬ 
periments of Raymond and Sessions (2007), Raymond and Fopez-Carrillo (2011, p. 147) 
summarized their findings as follows: "Environments cooler at low levels and warmer at 
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upper levels of order 1 K lower the elevation of maximum vertical mass flux from «10 
km to «5 km in their calculations, and by virtue of mass continuity, intensify the low-level 
inflow into the convection. In the presence of ambient vorticity, low-level convergence 
would contribute to the spinup of the system." In this paradigm for tropical cyclogenesis, 
the vorticity dynamics are slaved to the thermodynamics and mass continuity. 

This thesis focuses on the non-developing case of ex-Gaston (2010) and the developing case 
of Karl (2010) in the European Centre for Medium-Range Weather Forecasts (ECMWF) 
model. As discussed above, both of these systems were observed during the PREDICT 
experiment. Ex-Gaston is is arguably one of the most extensively observed non-developing 
tropical disturbances ever. During the PREDICT mission, there were five GV research 
flights over five days with a total of 109 dropsondes conducted and later ingested into the 
ECMWF model. 

Gjorgjievska and Raymond (2014, hereafter referred to as GR14) have proposed a different 
process for the failed redevelopment of ex-Gaston (discussed in detail Chapter 3.2). GR14 
question the structure of ex-Gaston’s pouch, as well as the role and efficacy of dry air in 
frustrating ex-Gaston’s redevelopment. This thesis answers two fundamental questions: 
(1) Did ex-Gaston have a robust pouch? (2) What was the role of dry-air, or anti-fuel, in 
ex-Gaston’s failed redevelopment? 

Unlike ex-Gaston, Karl underwent genesis and developed into a major hurricane. The 
PREDICT experiment began investigating this system nearly four days prior to genesis. 
This represents one of the earliest (with respect to time of genesis) and most observed 
developing systems. For Karl, there were six PREDICT research flights over five days with 
127 dropsondes conducted and later ingested into the EMCWF model. This thesis aims to 
establish the spin-up process for Karl. 

In closing this Introduction, it is useful to recall that a better understanding of the physics 
of tropical cyclogenesis would seem to be a necessary element for improving the prediction 
of storm development. An understanding of tropical cyclogenesis is of vital importance to 
the United States Navy. As humanitarian assistance and disaster relief becomes an increas¬ 
ingly important mission to promote regional stability, the navy will need to be effective at 
anticipating the vast destruction caused by tropical cyclones. A thorough understanding of 
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genesis will allow vital navy assets (e.g., carrier strike groups, amphibous ready groups, or 
the USNS Mercy or USNS Comfort) to anticipate when and where a tropical cyclone will 
form. This could provide a nearly 48-hour lead time for pre-positioning. In order to maxi¬ 
mize an efficient response while still protecting the fleet, decision makers rely on accurate 
and timely, environmental information. Navy Meteorology and Oceanography, therefore, 
needs the most current and innovative approaches for effective tasking, collection, process¬ 
ing, exploitation, and dissemination of tropical cyclone decision aids and recommendations 
to the fleet. 
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CHAPTER 2: 
Data Sources 


This study uses the European Centre for Medium-Range Weather Forecasts (ECMWF) 
analysis fields. The ECMWF model is a two-time-level, semi-implicit, semi-Lagrangian 
spectral model that solves the Euler equations of motions using the hydrostatic approxima¬ 
tion. The horizontal discretization uses a reduced Gaussian grid, while the vertical model is 
discretized using a finite-element scheme. The ECMWF model is able to resolve flow fea¬ 
tures to ~1°. The ECMWF model uses a four-dimensional variational analysis (4DVAR) 
system to incorporate observations into the model forecast. As summarized in the Intro¬ 
duction and in Montgomery et al. (2012), the PREDICT data consisting of the NSF GV 
dropsondes, together with the NASA DC-8 and NOAA G-IV and P3 dropsondes, were in¬ 
corporated into the ECMWF model. The resulting model analyses have a horizontal grid 
size of 0.25°, 25 vertical levels from 1 to 1000 hPa, and temporal output every 6-hours. 

Subgrid-scale process (e.g., radiation, clouds, and turbulent diffusion) are modeled through 
physical parameterization schemes of the model. ECMWF uses a moist convection scheme 
based on Tiedtke mass-flux approach (described in Tiedtke 1989) with a modified CAPE 
closure and an improved entrainment formulation. Bechtold et al. (2008) and Bechtold 
et al. (2014) describe the convection parameterization for ECMWF in detail. The radiation 
scheme uses the Rapid Radiation Transfer Model (Mlawer et al. 1997; Morcrette et al. 
2008) to compute the subgrid-scale radiative fluxes based on the predicted state of the 
atmosphere. ECMWF uses two different turbulent diffusions schemes: (1) a first order 
K-diffusion closure based on Monin-Obukhov similarity theory in the surface layer, and 
(2) above the surface layer an Eddy-Diffusivity Mass-Flux scheme is used for unstable 
boundary layers (Kohler et al. 2011). 

PREDICT mission scientists used model forecasts to design flight plans to investigate trop¬ 
ical disturbances during the experiment. The proposed flight plans attempted to perform 
dropsondes in a uniform grid with ~ 1° horizontal separation between drops (Figure 2.1). 
Because of flight restrictions from weather or international boundaries, some flight plans 
were altered and attaining a uniform distribution of dropsondes through the disturbance 
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was not possible (Montgomery et al. 2012). In the NSF supported PREDICT project, there 
were 5 research flights with 109 dropsondes conducted during ex-Gaston, and 6 research 
flights with 127 dropsondes conducted during Karl. Coincident with PREDICT, both the 
GRIP and IFEX experiments were also underway. The GRIP experiment used a NASA 
DC-8 aircraft to perform dropsondes into tropical disturbances. The IFEX experiment used 
G-IV and WP-3 aircrafts, and performed dropsondes in Karl and Matthew. The GRIP and 
IFEX data augmented the PREDICT dropsonde data set. The ECMWF model analyses 
used herein incorporate the PREDICT, the GRIP, and IFEX dropsonde data. 



(a) ex-Gaston (b) Karl 

Figure 2.1: Example flights paths for (a) ex-Gaston on 5 September and (b) Karl on 14 
September. Flight paths (red lines) are overlaid on minimum pressure contours (blue con¬ 
tours) derived from dropsonde data. Green lines indicate the length and direction of propa¬ 
gation of the sweet spot from the beginning of the research flight to the end of the research 
flight (from Smith and Montgomery 2012). 


For this thesis, the ECMWF analyses from 28 August to 11 September are used for ex- 
Gaston, and the analyses from 09 to 17 September are used for Karl. The sweet spot is 
determined through the following process: First, the wave-trough axis is determined in the 
ECMWF analyses as any line in an area of positive relative vorticity above 1 x 10 6 s 1 
along which meridional wind is zero. Second, the phase speed of the wave is resolved 
by determining the propagation of the meridional wind, total precipitable water, and/or 
relative humidity with Hovmoeller diagrams. Third, the sweet spot is then identified as the 
intersection of the wave-trough axis with the surface where zonal wind is equal to the phase 
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speed of the wave (critical latitude). 
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CHAPTER 3: 

Failed Redevelopment of Ex-Gaston 


Tropical Storm Gaston developed from an African easterly wave that moved westward from 
the African coast on 28 August 2010. The National Hurricane Center designated Gaston as 
a tropical storm at 1500 UTC 01 September. Despite being in a favorable environment with 
low shear and a SST of 28.5°C (Gjorgjievska and Raymond 2014), convection associated 
with Gaston quickly diminished, and NHC downgraded Gaston to a remnant low by 1800 
UTC 02 September. Convective activity increased on 03 September, however it did not 
re-organize and the system remained a remnant low (Figure 3.1). 



Figure 3.1: Track for ex-Gaston based on pouch center (black dots) as identified in the 6 
hour, ECMWF analysis data. Green squares show approximate times of PREDICT research 
flights over the disturbance. The red line indicates when the National Hurricane Center 
designated the disturbance as at least a tropical depression. The black-dashed line indicated 
when the disturbance was an incipient wave or remnant low. 


3.1 Ex-Gaston’s Environment 

This analysis begins by characterizing the vertical shear that affected Gaston’s pouch. 
“Deep-layer shear" and “pouch shear" are computed by taking the vector differential of 
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horizontal winds between the 200 hPa and 850 hPa levels, and between the 500 hPa and 
850 hPa levels, respectively, averaged over a a 3° x 3° box centered at the pouch center. 
The pouch-scale averaging is performed on a 3° x 3° box, centered on the circulation center 
as defined by the 700 hPa tracking level. 

For both the deep and pouch shear, the magnitude of the shear decreases rapidly from ~20 
ms -1 on 30 August to ~2 m s _1 on 02 September (Figure 3.2). During the same period, 
the direction of the deep and pouch shear shifts from easterly to westerly flow (Figure 
3.3). After 02 September, the magnitude of the shear (deep and pouch) increases to ~5 
ms -1 . The pouch shear direction slowly becomes more northerly by 05 September. The 
deep shear, though, rapidly changes direction from northeasterly to southwesterly from 
1200 UTC 02 September to 0000 UTC 03 September, in the ECMWF data. The deep 
shear returns to an easterly flow on 04 September. These shear results are consistent with 
the analysis of PREDICT data by Davis and Ahijevych (2012). The National Hurricane 
Center defines vertical shear of 12 m s' 1 as an upper limit for favorable conditions for 
tropical cyclogenesis. The magnitude of the vertical shear (typically 4-8 m s^ 1 ) for ex- 
Gaston, while below this heuristic limit for a SST of 28.5°C, does suggest lateral flow 
and a potential contribution of dipole-like distribution of vorticity from a non-advective 
flux (Haynes and McIntyre 1987; Raymond et al. 2014). This contribution could be a net 
increase or decrease of vorticity. 

The evolution of other pertinent variables is shown in Figure 3.4. In the subpanels of this 
figure, shown are time-height Hovmoeller diagrams of relative humidity, relative vorticity 
(£), and vertical mass flux at each level from averages taken over a 3° x 3° box (referred to 
as pouch scales). A similar analysis was done for a 1° x 1° box (referred to as sub-pouch 
scale) centered on the circulation center; trends were similar to those for the 3 degree box, 
but are not shown. This analysis uses the model vertical velocity in pressure coordinates, 
(U, to calculate the mass flux as pw — —G)/g, where p is density, w is vertical velocity in 
height coordinates, and g is the acceleration due to gravity. On both the pouch and sub¬ 
pouch scales in the ECMWF data, a layer of dry air above 600 hPa appears to penetrate the 
pouch region on 02 September, and that dry layer persists through the decline of the system 
(Figure 3.4). Coincident with the intrusion of the dry air are decreases in the system-scale 
relativity vorticity and mass flux. 
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Figure 3.2: Average of vertical wind shear magnitude is shown for the 200-850 hPa (red 
line) and the 500-850 hPa (blue line) levels. Shear decreased from ~ 20 m s _1 to ~ 2 m 1 
on 02 September 2010. Although the shear is below the nominal value of 12 m 1 , there 
is still persistent shear on the pouch. These results are consistent with previous studies by 
Davis and Ahijevych (2012). 


To gain insight into the apparent intrusion of dry air into the pouch beginning near 0600 
UTC 02 September (discussed above), this study first examines the flow topology of ex- 
Gaston using the dividing streamline methodology discussed in Riemer and Montgomery 
(2011). This methodology assumes for simplicity that the flow is steady in a translating 
frame. Although the observed flow will be shown later to have an important transient 
component, this technique can provide a first look into the existing flow topology around 
ex-Gaston’s pouch. Figure 3.5 shows the horizontal flow fields and calculated dividing 
streamlines at 1800 UTC 02 September from ECMWF analysis data at 500, 700, 850, and 
925 hPa levels when dry air was greatly impacting the pouch. 
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Figure 3.3: Direction of vertical shear analyzed over a 3° x 3° box centered on the moving 
sweet spot location. The vertical wind shear is shown for the 200-850 hPa (red line) and 
the 500-850 hPa (blue line) levels. The directions are compass directions and in a meteo¬ 
rological sense (e.g., 90° corresponds to easterly shear). The 500-850 hPa shear is mainly 
from the north on 02 September 2010. Analysis by RM12 shows that is a source region for 
dry air. 


At 700, 850, and 925 hPa, a hyperbolic point lies east of the circulation center. However, 
to the west of the circulation center, the pouch is open to the environment; thereby provid¬ 
ing a pathway for air parcels to enter the pouch (i.e., an “open pouch"). At 500 hPa the 
hyperbolic point is northwest of the circulation center, and the pouch is open to the east. 
This interpretation is consistent with the study of RM12 who performed a more compre¬ 
hensive study of the flow topology of ex-Gaston. In particular, RM12 analyzed Lagrangian 
coherent structures derived from particle trajectories, and found that lateral, dry air intru¬ 
sion occurred from 01 to 05 September (RM12 Figure 6 therein). These current findings 
at 700, 850, and 925 hPa are consistent with the detailed Lagrangian analysis of ex-Gaston 
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Figure 3.4: Time-height cross-section of system averaged quantities within a 3° x 3° box 
from 0000 UTC 30 August to 1800 UTC 06 September, (a) Relative humidity in percent, 
(b) Relative vorticity (£) with relative humidity contours of 40, 50, 60, and 70%. (c) 
Vertical mass flux with relative humidity contours of 40, 50, 60, and 70%. A dry layer near 
and above 600 hPa appears on 02 September and persists through 06 September. There are 
corresponding decreases in relative vorticity and mass flux after 02 September. 


by RM12. 

The structure of the pouch can be further studied by looking at the time-dependent nature 
of the flow by calculating hyperbolic trajectories (Samelson and Wiggins 2006). Hyper¬ 
bolic trajectories are trajectories of the time-dependent flow field that share the same linear 
stability properties as hyperbolic fixed points in time-independent flow. These hyperbolic 
trajectories have stable and unstable manifolds associated with them, and these manifolds 
control particle transport in time-dependent flow (Ide et al. 2002). 
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Figure 3.5: Dividing streamlines at 1800 UTC 02 September 2010. Dividing streamlines 
are shown for 500, 700, 850, and 925 hPa, and are overlaid on co-moving wind vectors at 
each level and total precipitable water. 


Figure 3.6 shows a time sequence of stable and unstable manifolds at 500 hPa (left column) 
and 700 hPa (right column) from 0000 UTC 01 September to 0000 UTC 03 September. 
Stable manifolds are indicated with red lines, and unstable manifolds are indicated with 
blue and cyan lines. For reference, a 3-degree radius circle around ex-Gaston’s diagnosed 
pouch center is indicated by the green circle. Throughout this 48-hour period, the stable 
manifold (red line) and an unstable manifold (blue line) intersect east of ex-Gaston’s pouch 
on the 500 and 700 hPa pressure surfaces. The intersection of these manifolds marks the 
location of a hyperbolic trajectory, and the persistence of these manifolds is indicative of 
the pouch having a barrier to intrusions from the northeast, east and southeast. In this case 
the manifolds comprise only part of a cat’s eye. At 700 hPa, the stable (red line) manifold 
also intersects the unstable (cyan line) manifold south of ex-Gaston’s pouch. This second 
intersection is not evident on the 500 hPa pressure level. At both the 500 hPa and 700 
hPa levels, there are no intersecting manifolds west of the pouch. No intersection implies 
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no additional hyperbolic trajectory, and leaves no way for boundaries to be topologically 
connected into a separatrix. 1 While there is a difference between the dividing streamline 
and Langrangian manifold analyses at 500 hPa, the two methods are consistent at the other 
pressure levels presented. The Langrangian manifold method is more complete and accu¬ 
rate by incorporating the time-dependent nature of the analyzed flow. 




0N 

50W 45W 40W 35W 30W 25W 20W 15W 







Figure 3.6: Lagrangian manifolds are overlaid on 9 e fields at 500 hPa (left column) and 
700 hPa (right column) from 01 to 03 September. Stable manifolds are in red, and un¬ 
stable manifolds are in blue and cyan. The abscissa is longitude and ordinate is latitude. 
The manifolds indicate that the pouch had a hyperbolic point to the east, but was open to 
environmental air to the west (figure provided by B. Rutherford, personal communication). 

*A separatrix is a flow partitioning boundary formed by connected segments of manifolds or material 
curves. 
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To identify the source region for the dry air that entered ex-Gaston’s pouch at 700, 600, 
500, and 400 hPa on 1800 UTC 02 September, backward trajectories were computed for 
particles seeded within a 3-degree radius of the pouch center. Trajectories were computed 
as in RM12 using a fourth-order Runge-Kutta method with a 15-minute intermediate time 
step and bi-cubic interpolation in both time and space on constant pressure surfaces. At 400 
hPa, particles that are within a 3-degree radius of the pouch at 1800 UTC 02 September 
originated to the north of the pouch (Figure 3.7). At the 500 hPa level, particles that are 
within a 3-degree radius of the pouch at 1800 UTC 02 September originated primarily 
northeast and southwest of the pouch (Figure 3.8). 



Figure 3.7: Particle trajectory locations at 400 hPa 1800 UTC 31 August are overlaid on 9 e 
(K). These trajectories are all within a radius of 3 degrees of the pouch center (green circle) 
by 1800 UTC, 02 September. The abscissa is longitude and ordinate is latitude (figure 
provided by B. Rutherford, personal communication). 

For the trajectories identified in the foregoing figures, it is of interest to document the 
evolution of pseudo-equivalent potential temperature, 9 e . For a moist air parcel, 0 e is ap¬ 
proximately materially conserved in the absence of mixing processes. On a given pressure 
surface, 9 e is a function of moisture and temperature and because of its tracer-like property 
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Figure 3.8: Particle trajectory locations at 500 hPa 1800 UTC 31 August are overlaid on 0 e 
(K). These trajectories are all within a radius of 3 degrees of the pouch center (green circle) 
by 1800 UTC 02 September. The abscissa is longitude and ordinate is latitude (figure 
provided by B. Rutherford, personal communication). 


and weak temperature gradient in the tropics, increases or decreases in 9 e along a constant 
pressure trajectory reflect primarily changes in moisture. For all calculations presented 
here, the 0 e definition as given by Bolton (1980) (his Eq. (43)) is used. 

Figure 3.9 summarizes the evolution of 0 e for the trajectories identified previously in Figure 
3.8 on the 500 hPa level. A colored point in the figure represents a snapshot of the particular 
particle’s 9 e and radial distance from the center of Gaston’s pouch. The colors range from 
brown to blue, with brown denoting the earliest time of 0000 UTC 31 August and blue 
denoting the latest time of 1800 UTC 2 September. The quasi-regular pattern of blue dots 
between 0 and 3 degrees radius is a manifestation of the initial seeding method for the 
backward trajectory calculation. 

Figure 3.9 shows that particles seeded within the nominal pouch radius of 3 degrees origi¬ 
nate from two distinct source regions (brown points) outside of the pouch. The two source 
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Figure 3.9: 9 e vs. radial distance from the center of ex-Gaston’s pouch. The colors range 
from brown to blue, with brown denoting the earliest time of 0000 UTC 31 August and 
blue denoting the latest time of 1800 UTC 02 September (figure provided by B. Rutherford, 
personal communication). 


regions are indicated by the red dots in Figure 3.8; one source region is located in an arch¬ 
like filament northeast of the pouch in a dry (low 0 e ~ 328 K) environment; the other is 
located west and southwest of the pouch in a relatively moist environment ( 6 e ~ 339 K). 
As these particles enter ex-Gaston’s pouch from 0000 UTC 31 August to 1800 UTC 02 
September, the moist particles remain relatively moist, and the dry particles remain rela¬ 
tively dry. The black slanted line in Figure 3.9 approximately differentiates these moist 
and dry trajectory paths, and its shallow slope indicates that dry air was not significantly 
moistened before entering the pouch. A similar analysis was performed for the 400 hPa 
level (not shown), and showed similar trends as the 500 hPa level. These results demon- 
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strate that dry air was entering Gaston’s pouch during this 66 hour period from 0000 UTC 
31 August to 1800 UTC 02 September. 


3.2 Implications of Dry Air and a Degraded Pouch 

The findings from the previous section showing dry air entering ex-Gaston’s pouch moti¬ 
vate an important question in its non-redevelopment: What was the role of dry air in the 
non-development? To understand the role of the dry air entering the pouch it is useful to 
review previous studies of the role of dry air on convection. A new hypothesis on the role 
of dry air in tropical cyclogenesis was articulated in the work of Smith and Montgomery 
(2012). The authors studied the convective environments of the tropical disturbances dur¬ 
ing the PREDICT experiment. They found that a prominent difference between developing 
and non-developing disturbances was the difference in 0 e between the surface and 3 km. 
Smith and Montgomery (2012) hypothesizd that entrained, dry air weakens the convective 
updraughts and thereby weakens the vortex-tube stretching of ambient and local cyclonic 
vorticity. Weakening of the convective updraughts implies a frustrated vorticity amplifi¬ 
cation process. The hypothesis of Smith and Montgomery (2012) stands in contrast to 
the traditional notion that dry air increases the strength of convective downdraughts and 
increases the low-level divergence that accompanies these downdraughts. 

In another study of convective environments, James and Markowski (2010) investigated 
the role of dry air aloft on deep convection. In their numerical study, they found in the 
low CAPE environments (1500 J kg -1 ) with a dry air layer of RH=70% near 700 hPa, 
the updraught mass flux was reduced throughout the depth of the troposphere, and the 
downdraught mass flux was either unchanged or reduced. 

In a related study, Kilroy and Smith (2012) (hereafter referred to as KS12) used a high- 
resolution cloud model to examine the role of dry air on rotating deep convection in tropical 
cyclogenesis. KS12 created an idealized sounding based on the ex-Gaston environment, 
and proceeded to modify the idealized sounding by injecting dry air into the mid-levels. 
They found through a series of experiments (summarized in their Table 2) that dry air 
aloft reduced the convective updraught strength and the vertical extent of the convective 
updraught. 

KS12 also used a “moist" and a “dry" profile from the PREDICT Experiment. The moist 
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profile was from 1820 UTC 05 September, and had a total precipitable water (TPW) of 65.2 
kg m 2 . The dry profile was from 1448 UTC 05 September, and TPW=43.5 kg m 2 . In the 
moist environment, KS12 found maximum convective updraught and downdraught veloc¬ 
ities of 34 and 10.9 m s 1 , respectively, and vertical extent above 10 km. However, in the 
dry environment the maximum updraught velocity and downdraught velocities were 11.4 
and 6.3 m s -1 , respectively, and the vertical extent was only ~ 7 km (Figure 3.10)). Dry air 
reduced both the updraught and the maximum vertical extent, while the downdraught ve¬ 
locities were only moderately reduced, consistent with findings from the experiments with 
the idealized soundings. These results showed that dry air reduces cloud buoyancy, thus 
making mass flux profiles weaker and shallower than in a moist environment, as well as 
making the convective updraught less effective in amplifying vertical vorticity over a deep 
tropospheric layer (Smith and Montgomery 2012). 
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Figure 3.10: Height-time series of maximum vertical velocities from the numerical experi¬ 
ments of KS12 using the (a) moist profile and the (b) dry profile from PREDICT. Red lines 
are positive (updraught) and blue-dashed lines are negative (downdraughts). Thin contours 
are 2 m s -1 and thick contours are 4 m s -1 . Thick black contours indicate values above 20 
m in5mE ! intervals (from KS12). 


The PREDICT data from 02 September are compared to the PREDICT data used by KS12 
(Figure 3.11). The moist profiles (Figures 3.11a and 3.1 lb) are from within the pouch, and 
show a moist environment in a layer from the surface to 400 hPa, and shallow, dry layer 
between ~300 - 400 hPa. From KS12, the moist sounding has CAPE = 3498 kg m 2 , CIN 
= 0 kg m~ 2 , and TPW = 65.2 kg m~ 2 (Smith and Montgomery 2012). A typical moist 
sounding from 02 September has CAPE = 1566 kg m~ 2 , CIN = 0 kg m~ 2 , and TPW = 59.9 
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kg m 2 (Smith and Montgomery 2012). 


The dry profiles (Figures 3.11c and 3.1 Id) were taken outside of ex-Gaston’s pouch. They 
both show a shallow moist layer below ~ 800 hPa, and dry through a deep layer above 
800 hPa. From KS12, the dry sounding has CAPE = 1145 kg m~ 2 , CIN = 21 kg m~ 2 , 
and TPW = 43.5 kg m 2 (Smith and Montgomery 2012). A typical dry sounding from 02 
September has CAPE = 525 kg m~ 2 , CIN = 75 kg m~ 2 , and TPW = 35.6 kg m~ 2 (Smith 
and Montgomery 2012). From the similarities in these profiles, one would expect similar 
results for 02 September as KS12 found for 05 September data (discussed above). 
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Figure 3.11: Comparison for moist (a) and dry (c) soundings from 05 September used by 
KS12 with typical moist (b) and dry (d) soundings from 02 September. Solid lines indicate 
temperature and dashed-lines indicate dew point temperature. Flight paths (white) for 05 
September (e) and 02 September (f) are overlaid on the RH field at 700 hPa (shading) and 
co-moving streamlines (black lines). Blue and red dots indicate the locations of the selected 
moist and dry soundings, respectively, for each day. 
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In total, the foregoing results suggest that the pouch was vulnerable to the environment with 
dry air penetrating the pouch and disrupting the amplification of vorticity. This analysis 
of the kinematic, dynamic, and thermodynamic structure of ex-Gaston in the ECMWF 
analysis, as well as the work of RM12 and Davis and Ahijevych (2012, Figure 9 therein), 
show that ex-Gaston’s pouch was misaligned, and vulnerable to environmental air as early 
as 02 September (Figure 3.4), the day of the first PREDICT research flight into this remnant 
low. This dry air results in divergence near the 600 hPa level, thus causing an expanding 
circulation loop at these levels. From Kelvin’s circulation theorem, as the loop expands 
the absolute vertical vorticity must decrease in order to conserve absolute circulation. A 
reduction in vorticity will create a compromised pouch, which will allow further intrusion 
of dry air and inhibit vorticity amplification. 

GR14 propose a different process that leads to the failed redevelopment of ex-Gaston. 
GR14 (p. 3076) hypothesize that the “severe decrease of the mid-level vortex observed 
between the period of Gaston 1 (02 September) and Gaston 2 (03 September) was a decid¬ 
ing factor for Gaston’s failure to re-intensify." It is important to note that GR14, RM12, and 
Davis and Ahijevych (2012) agree on the decay of the mid-level vortex. GR14 hypothesize 
that convection was suppressed by a strong trade wind inversion, and attribute the decrease 
in the mid-level vorticity to this strong trade inversion and corresponding structure of the 
vertical mass flux profile. GR14 argue that the trade wind inversion air causes the decrease 
in magnitude with height of the mass flux profile. However, this study contends that the 
intrusion of dry air at and above the 600 hPa level is responsible for the decrease in the 
mass flux profile. 

To address the hypothesized influence of the strong trade wind inversion, it proves useful to 
review the thermodynamic structure of the dropsonde data collected on 02 September. The 
PREDICT experiment released 19 dropsondes in ex-Gaston on this day. GR14 included 
drop numbers 2 through 14 (see Table 3.1) in their 3DVAR analysis and area-averaging 
schemes. Within their 4° x 4° analysis box (their Figure 8), only one drop (drop number 2, 
located northwest of ex-Gaston’s pouch) shows clear evidence of a temperature inversion 
(Figure 3.12), while 11 of the profiles show evidence of a dry layer above 600 hPa. It is 
unclear how one particular sounding could have such a hypothesized impact on the system- 
scale vorticity dynamics. 
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Table 3.1: Summary of Dropsondes from PREDICT Research Flight 9 (RF09) on 02 
September (after Smith and Montgomery 2012). 


Drop Num. 

Time 

(UTC) 

TPW 

(kgnw 2 ) 

CAPE 

(Jkg- 1 ) 

CIN 

(J kg” 1 ) 

1 

1532 

33.0 

478 

149 

2 

1544 

48.2 

196 

95 

3 

1555 

53.6 

24 

142 

4 

1605 

61.1 

688 

6 

5 

1614 

62.8 

706 

0 

6 

1624 

57.5 

1047 

9 

7 

1637 

58.4 

612 

29 

8 

1647 

63.0 

1707 

0 

9 

1654 

65.9 

654 

11 

10 

1703 

67.1 

1649 

0 

11 

1713 

59.9 

1566 

0 

12 

1723 

57.5 

605 

5 

13 

1733 

56.7 

2 

158 

14 

1745 

55.3 

0 

328 

15 

1755 

53.8 

114 

110 

16 

1808 

51.1 

1155 

14 

17 

1818 

35.6 

525 

75 

18 

1830 

36.1 

285 

143 

19 

1843 

38.2 

101 

155 


A study of GR14 suggests that these authors appeared to overlook the implications of Davis 
and Ahijevych (2012) findings of a vertically sheared pouch and RM12’s findings of dry 
air mixing into ex-Gaston’s pouch between 01 and 03 September (RM12’s Figure 6), a 
time period spanning the first day of PREDICT observations (02 September). While GR14 
acknowledge the role of a transient flow component in causing a reversal in the sign of 
the vorticity tendency, they do not recognize the role of this component in modulating the 
transport of dry air into ex-Gaston. Thus, GR14 appear also to misinterpret the results of 
Smith and Montgomery (2012) and RM12: GR14 (p. 3076-3077) imply that ex-Gaston’s 
pouch was “robust" (i.e., closed) on 02 September, and therefore unlikely that “dry air 
might have been drawn into the core of Gaston." 

The offered hypothesis of Gaston’s non-redevelopment described above stands somewhat 
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Figure 3.12: PREDICT flight path and dropsonde locations (white line with white points) 
overlaid on 700 hPa relative humidity (shading) and co-moving streamlines from ECMWF 
analysis data at 1800 UTC 02 September. The gray Xs indicate locations where dropsonde 
data show evidence of an inversion. The black box corresponds to GR14 Figure 8. The 
expanded box shows the skew T - In P diagrams from the 02 September PREDICT flight 
dropsondes within the GR14 area of interest. Only the profile in the northwest comer in 
the black box shows evidence of an inversion. 


in contrast to the alternative hypothesis by Gjorgjievska and Raymond (2014). The data 
shown herein supports the view that dry air penetrated the pouch before the first flight into 
ex-Gaston and disrupted the amplification of vorticity at those levels where dry air intruded. 
GR14 agree with this hypothesis of a compromised pouch and dry air intrusion, but only 
after 04 September. 

The data from PREDICT research flight 9 on 02 September does show evidence for an in¬ 
version outside the GR14 analysis box for ex-Gaston’s pouch (Figure 3.12). To investigate 
the influence of the dry, trade inversion air west of ex-Gaston’s pouch on 02 September, 
a forward trajectory analysis (Figure 3.12, grey Xs west of the black box) was performed. 
Particles were seeded on the 850 hPa pressure level west of the sweet spot (where the 
PREDICT data show a temperature inversion) at 1200 UTC 02 September and integrated 
forward to 0000 UTC 04 September. Nearly all of the particles are located outside of a 3 
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degree radius of the pouch center by 0000 UTC 04 September. This analysis shows that the 
observed dry, trade inversion air on 02 September does not enter the pouch, and has little 
influence on the non-development of ex-Gaston. 
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Figure 3.13: Particle trajectory locations at 850 hPa at 0000 UTC 04 September overlaid 
on 9 e (K). Particles were seeded west of ex-Gaston’s pouch at 1200 UTC 02 September, 
within a suspected trade wind inversion layer. These trajectories are all outside of a 3 
degree radius of the pouch center (green circle) on 04 September. The abscissa is longitude 
and ordinate is latitude (figure provided by B. Rutherford, personal communication). 
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CHAPTER 4: 

The Development of Hurricane Karl 


Hurricane Karl developed within a region of the intertropical convergence zone (ITCZ) 
near the northern coast of South America. South westerlies off the continent helped to en¬ 
hance the ITCZ and increase convection between 50° and 60° W on 08 September (Mont¬ 
gomery et al. 2012). An easterly wave formed near this enhanced region of the ITCZ and 
continued to propagate westwards into the Caribbean. The tropical low strengthened to 
a tropical storm on 2100 UTC 14 September, and made landfall along Mexico’s Yucatan 
coast. Karl then strengthened to a category 3 hurricane over the Bay of Campeche on 17 
September, and made landfall near Vera Cruz, Mexico (Figure 4.1). Although there was 



Figure 4.1: Track for Karl based on pouch center (black dots) as identified in the 6 hour, 
ECMWF analysis data. Green squares show approximate times of PREDICT research 
flights over the disturbance. The red line indicates when the National Hurricane Center 
designated the disturbance as at least a tropical depression. The black-dashed line indicated 
when the disturbance was an incipient wave or remnant low. 

recurrent deep, mesoscale convection in and around the pre-Karl disturbance, it was not 
until approximately 1200 UTC 14 September that there was organized and persistent deep 
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convection to support the transition of the pre-Karl wave-pouch disturbance into a tropi¬ 
cal cyclone. Figure 4.2 summarizes the evolution of the convection associated with Karl’s 
development. 



Figure 4.2: 85 GHz montage for Karl’s track from 10 to 17 September. Warmer colors 
indicate active convection (from Montgomery et al. 2012). 


Despite favorable conditions for tropical cyclone formation, Karl’s development was not 
a forgone conclusion. PREDICT mission scientists observed periods of deep convection 
that nearly coincided with the diurnal cycle (Davis 2010a,b). Davis and Ahijevych (2012) 
presented a time-radius diagram of cloud-top temperatures derived from Geostationary Op¬ 
erational Environmental Satellite (GOES) IR data during Karl’s development (Figure 4.3). 
This data indicates that cold cloud signature decreases simultaneously at all radii. Davis 
and Ahijevych (2012) notes that the convection (or cold cloud signature) initiates between 
1200 and 1500 UTC. These times correspond to the diurnal peak in convection for the 
Caribbean (Hirose et al. 2008; Sorooshian et al. 2002). After genesis, Figure 4.3 shows 
there is a persistent cold cloud signature extending to a radius of ~150 km (Davis and 
Ahijevych 2012). 
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Figure 4.3: Time-radius diagram of cloud-top temperatures derived from GOES IR data. 
Temperatures are the 75th percentile temperature within each radial ring of 20 km. The 
dashed line indicates the time of genesis. The black arrows indicate the times of GV PRE¬ 
DICT flights and the red arrows indicate times of research flights by NOAA or NASA (from 
Davis and Ahijevych 2012). 

4.1 Pre-Karl’s Environment 

This analysis of Karl’s genesis begins by characterizing its development environment. Sim¬ 
ilar to the ex-Gaston analysis, the “Deep-layer shear" and “pouch shear" are computed by 
taking the vector differential of horizontal winds between the 200 hPa and 850 hPa levels, 
and between the 500 hPa and 850 hPa levels, respectively, averaged over a 3° x 3° box 
centered at the pouch center. The pouch-scale averaging is performed on a 3° x 3° box, 
centered on the circulation center as defined by the 700 hPa tracking level. 

The pouch shear remains low (~2-8 ms -1 ) from 09 to 17 September. The deep-layer shear 
is similar to the pouch scale, except for a period from 0000 UTC 10 September to 1200 
UTC 11 September when the magnitude of vertical shear increases to ~13 m s” 1 (Figure 
4.4). For both the deep-layer and pouch shear, the direction is predominantly from the east 
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and southeast (Figure 4.5). From 1800 UTC 11 September to 1800 UTC 12 September, the 
direction of the pouch shear rotates clockwise from the east to the west, and then back to 
the east. These shear results are consistent with the analysis of PREDICT data by Davis 
and Ahijevych (2012, Figure 10 therein). Karl eventually developed in a favorable, low 
shear environment. 



Figure 4.4: Average of vertical wind shear magnitude is shown for the 200-850 hPa (red 
line) and the 500-850 hPa (blue line) levels. 
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Figure 4.5: Direction of vertical shear analyzed over a 3° x 3° box centered on the moving 
sweet spot location. The vertical wind shear is shown for the 200-850 hPa (red line) and the 
500-850 hPa (blue line) levels. The directions are compass directions in a meteorological 
sense (e.g., 90° corresponds to easterly shear). 

The evolution of relative humidity, relative vorticity (£), and vertical mass flux is analyzed 
using time-height Hovmoeller diagrams (Figure 4.6). The quantities were averaged over 
3° x 3° boxes centered on the sweet spot as defined by the 700 hPa tracking level. As 
in Chapter 3, this analysis uses the model vertical velocity in pressure coordinates, (0 , to 
calculate the mass flux as pw — —t'o/g , where p is density, w is vertical velocity in height 
coordinates, and g is the acceleration due to gravity. The relative humidity shows that 
while there was a dry layer above 600 hPa from 0000 UTC 09 September to 0000 UTC 11 
September, pre-Karl was in a moist environment with average relative humidities generally 
exceeding 80% at all levels. Calculations of 0 e (as defined by Bolton 1980, his Eq. 43) 
averaged over a 3° x 3° box are a complimentary way to demonstrate that Karl became 
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more moist throughout the mid-levels from 09 to 14 September (Figure 4.7). 





Figure 4.6: Time-height cross-section of system averaged quantities within a 3° x 3° box 
from 0000 UTC 09 September to 0000 UTC 17 September, (a) Relativity humidity in 
percent, (b) Relative vorticity (£) with relative humidity contours of 40, 50, 60, and 70%. 
(c) Vertical mass flux with relative humidity contours of 40, 50, 60, and 70%. 
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Figure 4.7: Time-height plot of 0 e averaged over a 3° x 3° box from 09 to 17 September. 
Higher values of 0 e indicate greater moisture. 

The relative vorticity between 800 and 1000 hPa at 1200 UTC 09 September (Figure 4.6 
shows that pre-Karl started to spin-up at these levels. The spin-up extends to the mid¬ 
levels (near ~500 hPa) by 0000 UTC 10 September. The pre-Karl disturbance then under¬ 
goes episodic spin-up and spin-down events before undergoing genesis on 1800 UTC 14 
September. At low-levels (850 and 950 hPa), the system spins-down near 0000 UTC 12 
and 0000 UTC 13 September (Figure 4.8). The 700 hPa £ is steady from 1200 UTC 10 
September to 0000 UTC 13 September. By 1800 UTC 14 September, Karl is spinning-up 
at all levels. 


37 










Figure 4.8: Evolution of £ averaged over 3° x 3° box at the 500, 700, 850, and 950 hPa 
pressure levels. The time is from 0000 UTC 09 September to 0000 UTC 17 September. 
Genesis occurred at approximately 1800 UTC 14 September. 


In order to investigate the aggregation and amplification of relative vorticity, the distribution 
of £ was calculated as a function of radius (Figure 4.9). Initially, the £ was averaged over 
a 1° x 1° box centered on the tracking level sweet spot (700 hPa). £ was then calculated in 
concentric square annuli with a 0.5° widths for total box sizes of 2° x 2° to 4° x 4°. Wang 
(2012) noted that averaging £ over pouch scale areas will tend to reduce the calculated 
magnitude of £. The averaging scheme used in this calculation is sub-pouch scale and 
should adequately represent the characteristic distribution of £ as a function of radius. 
At both the 700 hPa (Figure 4.9a) and 950 hPa levels (Figure 4.9b) there is a moderate 
amplification of vorticity from 0000 UTC 11 September to 1200 UTC 11 September. At 
700 hPa, the pulsing of £ is clear within a box size of 2° x 2°. The pulsing is much less 


38 








pronounced at the 950 hPa level. This could be a result of the averaging algorithm using 
the 700 hPa sweet spot position to define the box locations for all levels. As discussed 
by Davis and Ahijevych (2012) and later in this thesis, pre-Karl’s pouch was vertically 
misaligned until late 13 September. Therefore, the box may not be centered on the sweet 
spot for levels other than the tracking level. However, at 1800 UTC 14 September, there is 
a rapid amplification of £. As discussed above, there was persistent, deep convection near 
pre-Karl’s sweet spot at 1200 UTC 14 September, and pre-Karl underwent a transition to a 
tropical storm by 1800 UTC 14 September. 



Box size centered on sweet spot (deg) Box size centered on sweet spot (deg) 


(a) 700 hPa (b) 950 hPa 

Figure 4.9: Time-radius plots of average £. £ was averaged over a 1° x 1° box centered 
on the 700 hPa (tracking level) sweet spot, and subsequently averaged within 0.5° wide 
concentric square annuli. Results are shown for the 700 hPa (a) and 950 hPa (b) levels. 


The evolution of mass flux (Figure 4.6) shows a pulsing of upward mass flux (warm colors) 
until the system fully develops into a tropical cyclone. The pulses of upward mass flux 
are indicative of the intermittent deep convection observed during the PREDICT research 
flights into this disturbance (Montgomery et al. 2012). Mature mesoscale convective sys¬ 
tems (MCS) west of the 900 hPa circulation center (Davis and Ahijevych 2012) support the 
strong ascent observed in the ECMWF analysis data above 600 hPa at 1200 UTC 11 and 
12 September. Sustained, deep convective activity associated with Karl begins at approxi¬ 
mately 1200 UTC 14 September. 
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4.2 Circulation and Vorticity Dynamics 

Following the treatments of Haynes and McIntyre (1987), Raymond and Lopez-Carrillo 
(2011), and Lussier III et al. (2014), the equation describing the evolution of vertical vor¬ 
ticity can be written as follows: 


dr) 

~di 


—V h • Z + k ■ V/,0 x V/,II, 


(4.1) 


where 77 is the absolute vertical vorticity (/ + £), the h subscript denotes gradients on 

horizontal surfaces, 0 is the potential temperature, and FI is the Exner function. The Exner 

'k 

function is defined as II = , where p is the pressure, p () is a reference pressure, Rj 

is the gas constant for dry air, and c p is the heat capacity of dry air at a constant pressure. 
The vorticity flux, Z, is defined as (Lussier III et al. 2014): 


Z - Zi + Z 2 + Z 3 - r) z v h - T) h v z + k x F, (4.2) 


where a subscript of z denotes vertical vector components, v/ 7 are horizontal components 
of velocity, and F represents horizontal friction and other sub grid-scale forcing. The term 
k ■ V /; 0 x V/jII is the vertical component of the baroclinic generation of vorticity, and is 
generally neglected because of its small magnitude in the tropics (Raymond and Lopez- 
Carrillo 2011). In order to obtain an expression for the circulation tendency, one integrates 
Equation (4.1) (neglecting the relatively small baroclinic generation term) over a horizontal 
surface area, and uses Gauss’s divergence theorem to write: 


dV a 

dt 


v n T) z dl + (p r) n v z dl + j> F t dl , 


(4.3) 


where T u is the absolute circualtion, n denotes the outward normal component of a vector, 
t is the tangential component of a vector, and dl is line element along a circuit that defines 
the perimeter of the integration area. 


As described in Raymond et al. (2014), the first term of the right hand side of Equation (4.3) 
represents the net flux of vertical vorticity into the closed loop. The absolute circulation 
will not change unless there is a convergence of positive vertical vorticity (Figure 4.10, 
left panel). When this happens, the absolute circulation of the loop (and the area-averaged 
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vertical vorticity within the loop) will increase with time (Haynes and McIntyre 1987). 
The second term is a non-advective flux of vertical vorticity, and is commonly referred to 
as the vortex tilting-like term. This term can change the sign of vorticity (Raymond and 
Lopez-Carrillo 2011). Upward vertical motion in the region of a horizontal flux of vertical 
vorticity can distort this surface of equal vorticity upward (Figure 4.10, center panel). This 
mechanism creates a region of positive vertical vorticity and a region of negative vertical 
vorticity, and results in a vorticity tendency dipole (Raymond et al. 2014). The third term 
represents friction (Figure 4.10, right panel) or, for this analysis, any other subgrid-scale 
forcing in the ECMWF analysis data. 



Figure 4.10: Depiction of vorticity fluxes (from Raymond et al. 2014). 


In order to understand the circulation tendency and corresponding circulation budget for 
pre-Karl, the circulation and the advective and non-advective fluxes were calculated on 
each ECMWF pressure level above and including 950 hPa. Each quantity was calculated 
and averaged over boxes from 0.5° to 8.0° in 0.5° increments at each 6-hour time step of 
the ECMWF analysis data. The absolute circulation tendency was calculated by subtracting 
the previous 6-hour circulation field, from the current circulation field. As an example, the 
tendency at 1800 UTC 12 September is obtained by subtracting 1200 UTC from 1800 UTC 
12 September. 

During the development of Karl, there were episodes of spin-up and spin-down at the mid¬ 
dle and low-levels (Figures 4.6 and 4.8). As an example, the period from 1800 UTC 11 
September to 0000 UTC 12 September shows negative circulation tendency and conver¬ 
gence of absolute vorticity tendency below 700 hPa, but positive tendencies above 700 hPa 
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(Figure 4.11). From 0000 UTC 12 September to 0600 UTC 12 September, the tendencies 
are opposite with positive circulation and convergence of absolute vorticity tendencies be¬ 
low 600 hPa, but negative tendencies above 600 hPa (Figure 4.12). The contribution from 
the non-advective (tilting-like) tendency is generally found here to be minimal except in a 
shallow layer at large box sizes. 


Circulation Tendency 



Box size centered on sweet spot (deg) 

Figure 4.11: Time-height cross-sections of circulation (top), absolute vorticity convergence 
(middle), and vortex tilting (bottom) tendencies (shading) from 1800 UTC 11 September 
to 0000 UTC 12 September. The dashed line represents the zero contour level on each 
panel. This time interval is chosen to illustrate a low-level spin-down phase as indicated in 
Figures 4.6 and 4.8. 
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Circulation Tendency 



1 2 3 4 5 6 7 8 

Box size centered on sweet spot (deg) 


Figure 4.12: Time-height cross-sections of circulation (top), absolute vorticity convergence 
(middle), and vortex tilting (bottom) tendencies (shading) from 0000 UTC to 0600 UTC 
12 September. The dashed line represents the zero contour level on each panel. This time 
interval is chosen to illustrate a low-level spin-up phase as indicated in Figures 4.6 and 4.8. 

The pulsing of the circulation tendency continues until approximately 1800 UTC 14 
September. At this time, the circulation and convergence of absolute vorticity tendencies 
are positive through a deep layer (Figure 4.13). This positive trend continues until Karl 
develops into a category 3 hurricane. 
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Figure 4.13: Time-height cross-sections of circulation (top), absolute vorticity convergence 
(middle), and vortex tilting (bottom) tendencies (shading) from 1200 UTC to 1800 UTC 
14 September. The dashed line represents the zero contour level on each panel. NHC 
designated Karl a tropical storm at 2100 UTC 14 September. 

At the time of the PREDICT mission, the science director provided two hypotheses for 
the pulsing nature of the relative vorticity and circulation tendency (Montgomery 2010b): 
(1) the diurnal cycle for convective activity or (2) a boundary layer recovery time-scale 
related to convective activity, itself. The second hypothesis assumes that each convec¬ 
tive maximum is followed by a period of convective and mesoscale downdraughts. These 
downdraughts tend to flood the boundary layer with low 0 e air. A period of boundary layer 
recovery is needed to increase 0 e in the boundary layer back to pre-low 0 e anomaly levels. 
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In the absence of such convective updraughts and downdraughts during this recovery pe¬ 
riod, Montgomery (2010b) suggests that a simple slab boundary layer model can be used 
to estimate the boundary layer recovery time-scale as: 


T r = 


h 


(4.4) 


where T r is the recovery time, h is the depth of the boundary layer, Ck is the enthalpy 
exchange coefficient, and v is the near-surface wind speed. Using values of Ck — 1.2 x 10~ 3 
(Bell 2010), h = 0.5 km, and v = 8 ms -1 , one obtains T r «U5 h. Although this is an 
oversimplified model, Montgomery (2010b) finds this recovery time-scale is comparable to 
the period of observed convective activity and pulses in the relative vorticity and circulation. 
Montgomery (2010b) concludes that the periodic nature of the deep convective activity is 
result a of both the diurnal cycle in the Caribbean, as well as the boundary layer recovery 
time-scale. 


4.3 Karl’s Spin-Up 

The circulation tendency results show an overall trend of Karl spinning-up from 10 to 17 
September. A concomitant analysis is to study the evolution of the azimuthally-averaged 
absolute angular momentum surfaces and azimuthally averaged tangential velocity. The 
absolute angular momentum per unit mass, M is defined as follows: 

M = rv+^/r 2 , (4.5) 

where r is the radial distance to the sweet spot, v is the azimuthally averaged tangential 
velocity, and / is the Coriolis parameter at the sweet spot. In the absence of any torques 
(e.g., from eddy forcing, friction, or other subgrid-scale forcing), the absolute angular mo¬ 
mentum is materially conserved (i.e., DM/Dt = 0). Equation (4.5) can also be re-written 
in terms of the azimuthally averaged tangential velocity as: 

M 1 

v =- -fr. (4.6) 

r 2 

When M is materially conserved, v is directly related to the radial distance, r, of the M 
surface from the sweet spot. As an air parcel gets closer to the sweet spot, the azimuthally 
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averaged tangential velocity must also increase. Therefore, a necessary condition for trop¬ 
ical cyclone spin-up, in an inviscid and axisymmetric framework, is radial inflow toward 
the sweet spot (Montgomery and Smith 2014). 

Absolute angular momentum is approximately materially conserved above the frictional 
boundary layer. Deep convection near the sweet spot induces inflow above the boundary 
layer. Montgomery and Smith (2014) found in a number of idealized numerical simula¬ 
tions that above the frictional boundary layer, although the inflow is weak, the inflow is 
persistent. This persistent inflow draws M surfaces inwards, thus increasing the tangential 
winds above the boundary layer (Figure 4.14a,b). 

Within the frictional boundary layer, M is not materially conserved. Friction reduces the 
tangential velocity in this layer. While the tangential velocity is decreasing, the radial com¬ 
ponent of pressure gradient remains unchanged, and results in an inward directed agradient 
force. This agradient force generates a strong radial inflow (Figue 4.14c) in the vortex 
boundary layer. Because of the frictional effects and the strong radial inflow, the M sur¬ 
faces will tend to have a vertical gradient near the surface (Montgomery and Smith 2014). 
Figure 4.14d provides a summary depiction of this spin-up process involving both the deep 
convection and surface friction. 
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M reduced by friction, but strong convergence means 
that air parcels can reach small r quickly giving large v 


Figure 4.14: Absolute angular momentum surfaces (red and pink contours) at (a) initial 
time, and (b) 48 hours later from an idealized numerical experiment. The cyan contour 
indicates an M value of 20 x 10 5 m 2 s _1 , and the yellow region indicates values greater than 
32 x 10 5 m 2 s . Black contours in (b) are tangential velocity, and blue contours in (c) are 
radial velocities at 48 hours. The black contour in (c) delineates regions where conditions 
are met for axisymmetric instability (Montgomery and Smith 2014). (d) Depiction of the 
spin-up process (from Montgomery and Smith 2014). 


To study the evolution of M surfaces (calculated as M = T a /2n) for Karl, vertical profiles 
of M and v were made at each 6-hour timestep from 0000 UTC 10 September to 0000 UTC 
17 September. These calculations assumed a constant / at each time interval and corre¬ 
sponded with the latitude of the tracking level sweet spot. The general trend shows the M 
surfaces moving inward throughout Karl’s development (Figure 4.16). Initially, the radial 
inflow of air draws M surfaces in faster at larger radii. By 1800 UTC 13 September, the M 
surfaces are vertically uniform through a layer between 400 to 950 hPa. After 1800 UTC 
13 September, the contours are continually drawn in through a deep layer. Within the fric¬ 
tional boundary layer however, the combination of friction and strong radial inflow creates 
the vertical gradient in the M surfaces calculated from the ECMWF analysis data. This 
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spin-up process is consistent with the intensification paradigm described in Montgomery 
and Smith (2014) and summarized conceptually in Figure 4.14d. 

In a study of the structural evolution of Karl using the PREDICT data, Davis and Ahi- 
jevych (2012) observed a misalignment between the midlevel and low-level circulation 
centers (Figure 4.15). By 14 September, the circulation centers become aligned and fur¬ 
ther enhance a favorable formation environment as described in DMW09. With a favorable 
environment and pouch, pre-Karl transitioned into tropical storm Karl by 2100 UTC 14 
September. 



-90 -85 -80 -75 -70 -65 -60 -55 


Longitude (W) O 900 hPa O 500 hPa 

Figure 4.15: Pre-Karl circulation centers at 500 hPa (blue) and 900 hPa (red). The date and 
time of the analysis is indicated by the blue labels. The size of the circle is related to the 
spatial extent of the circulation (from Davis and Ahijevych 2012). 


Combining the ideas of the inward-moving M surfaces and the approximate material con¬ 
servation of M, Equation (4.6) shows that the tangential winds intensify. This amplification 
is, in fact, observed for Karl (Figure 4.16). Early in Karl’s development, the maximum 
tangential wind is ~11 m s -1 in a layer between 700 and 950 hPa. As the system develops, 
the tangential wind increases at all radii through a uniform, deep layer from 400 to 950 hPa. 
In general, the evolution of M surfaces and azimuthally averaged tangential wind implies 
that Karl underwent a low-level spin-up. This azimuthally-averaged perspective of Karl’s 
development is similar to the findings of Lussier III et al. (2014) for Typhoon Nuri, and 
consistent with the intensification paradigm presented by Montgomery and Smith (2014). 
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Figure 4.16: Evolution of average tangential velocity (shading) and angular momentum 
surfaces (contours) from 1800 UTC 11 September to 1800 UTC 16 September in 24 hour 
intervals as examples. The ordinate is pressure from 200 to 950 hPa and the abscissa is the 
size of a box centered on the sweet spot. Contours are shown for 0.5 x 10 6 m 2 s 'intervals 
with the 1, 3, 5, and 7x 10 6 m 2 contours highlighted in white. 
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CHAPTER 5: 
Conclusions 


The marsupial paradigm for tropical cyclogenesis hypothesizes that convective activity is 
favored within the protective region of a robust pouch (DMW09). Barring unfavorable 
circumstances like dry air intrusion or persistent vertical misalignment of the pouch, the 
theory hypothesizes that a developing system will exhibit a progressive moistening of the 
pouch by convection. Furthermore, vortical convective activity and convergence of vortic- 
ity near the sweet spot will amplify and aggregate vorticity in the pouch. Amplification and 
aggregation of vorticity will lead to the genesis of a tropical cyclone. 

In this thesis, evidence is presented to suggest that both the non-developing case of ex- 
Gaston and developing case of Hurricane Karl support the marsupial paradigm for tropical 
cyclogensis. Using the ECMWF analysis data, this thesis demonstrates that ex-Gaston did 
not have a robust pouch and was open to the intrusion of environmental air at the mid to 
upper levels on all days of the aircraft observations of this disturbance during PREDICT. 
Lagrangian trajectory and manifold calculations using ECMWF analyses show that dry air 
did indeed penetrate the pouch. These findings support the hypothesis that entrained, dry 
air near 600 hPa inhibited convective updraughts and vertical mass flux, which in turn leads 
to a reduction in vorticity and a compromised pouch at these middle levels. A compromised 
pouch allows further intrusion of dry air and inhibits subsequent vorticity amplification, as 
articulated in the work of Smith and Montgomery (2012). The findings presented herein 
support the original hypothesis that ex-Gaston’s degraded pouch further led to the non¬ 
redevelopment of the system by limiting the amplification of vorticity and not providing 
a protected environment for sufficient vorticity aggregation, consistent with the marsupial 
paradigm of tropical cyclogenesis described by DMW09. 

Unlike ex-Gaston, pre-Karl was embedded within a moist environment and developed into 
a tropical storm. Although Karl was not an African easterly wave (like ex-Gaston), a wave 
interaction with a mesoscale disturbance within the ITCZ led to the formation of a pouch 
(RM12). Despite forming in a favorable environment (warm SST, moist air, and low ver¬ 
tical shear), Karl had events of punctuated convection that resulted in both negative and 
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positive tendencies in the circulation throughout its development until its eventual genesis 
on 14 September. Davis and Ahijevych (2012) described the initial, vertical misalignment 
of Karl’s pouch that helped to impede Karl’s development. As this current study shows, by 
1800 UTC 14 September (coincident with the vertical alignment of the pouch), convection 
within Karl’s pouch was able to amplify relative vorticity; leading to a positive circulation 
tendency, an inward flux of angular momentum surfaces, and enhanced tangential wind. 
Karl’s development is consistent with the three hypotheses of tropical cyclogenesis posited 
by DMW09 and summarized in Chapter 1 of this thesis, as well as with the intensification 
model described by Montgomery and Smith (2014). 
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